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Abstract: we experimentally demonstrate an all-fiber spectrometer 
consisting of a photonic bandgap fiber bundle and a black and white 
CCD camera. Photonic crystal fibers used in this work are the large sohd 
core all-plastic Bragg fibers designed for operation in the visible spectral 
range and featuring bandgaps of 60nm - 180nm-wide. 100 Bragg fibers 
were chosen to have complimentary and partially overlapping bandgaps 
covering a 400nm-840nm spectral range. The fiber bundle used in our 
work is equivalent in its function to a set of 100 optical filters densely 
packed in the area of ~lcm^. Black and white CCD camera is then used 
to capture spectrally "binned" image of the incoming hght at the output 
facet of a fiber bundle. To reconstruct the test spectrum from a single 
CCD image we developed an algorithm based on pseudo-inversion of 
the spectrometer transmission matrix. We then study resolution limit of 
this spectroscopic system by testing its performance using spectrally 
narrow test peaks (FWHM 5nm-25nm) centered at various positions 
within the 450nm-700nm spectral interval. We find that the peak center 
wavelength can always be reconstructed within several percent of its true 
value regardless of the peak width or position. Moreover, we 
demonstrate that although the widths of the individual Bragg fiber 
bandgaps are quite large (>60nm) the spectroscopic system has a 
resolution hmit of ~30nm. We conclude by showing theoretically, that, 
in principle, presented spectroscopic system can resolve much narrower 
peaks down to several nm in width; this would, however, require a 
stringent control of the experimental errors during measurement and 
calibration. We believe that photonic bandgap fiber bundle-based 
spectrometers have a potential to become an important technology in 
multispectral imaging because of their simplicity (lack of moving parts), 
instantaneous response, and high degree of integration. 
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1. Introduction 

Development of the multispectral imaging systems presents an active research field in rapid 
development. Inherent advantages of such systems are in the real-time monitoring, 
cost-effectiveness and high accuracy. Imaging spectrometers have been demonstrated in many 
applications including chemical component analysis [1], remote sensing [1,2], astronomy [3], 
and biochemistry [4]. Moreover, multispectral imaging systems have been successfully 
demonstrated both in the visible and the near inferred spectral ranges [5-9]. 

A typical multispectral imaging system includes a high sensitivity, high signal-to-noise 
ratio photo-electric conversion device, such as a CCD array; optical filters or dispersive 
elements which are either used to provide the reference spectra or to filter-out a specific 
spectral component; as well as a broadband illuminant (light source). In a typical imaging 
spectrometer, spectral information is acquired indirectly by taking several images in the 
complimentary, while possibly overlapping, spectral ranges. A set of intensity images is then 
converted into spectra by utilizing various interpretation algorithms such as pseudo-inverse 
estimation (PSE), principal-component analysis (PGA), and nonlinear fitting [5,6]. 

Recently, several optical fiber-based imaging spectrometers have been reported [7,8]. In 
[7], for example, the authors used at the input end a 10x10 bundle of the mid infrared fibers 
which were then separated into a linear 1 x 100 array at the output end. The output beams were 
then simultaneously dispersed using a single Bragg grating and imaged on a CCD array. This 
system, thus, functions as a 10x10 spectral imaging array. In another implementation [8], a 
multiple range spectrophotometer was implemented using an optical fiber probe. In that 
system, linear displacement of the light source was projected onto a fiber array (bundle) so 
that only a single fiber was lit for a given position of the light source. The fiber bundle was 
then interrogated using a spectrometer, thus enabling spatially resolved real-time data 
acquisition. Note that in both examples, optical fibers were used exclusively as broadband 
optical guides, while a separate spectrometer was still required for the spectra acquisition. 

The goal of this paper is to demonstrate that the function of a spectrometer can be 
implemented directly inside a single fiber bimdle. Then a complete imaging spectrometer can 
be realised using an array of such fiber bundles interrogated by a linear or a 2D CCD sensor. 
Such an approach would allow getting rid of costly and slow traditional spectrometers based 
on moving gratings, thus leading to significant cost savings and increased acquisition speed. 
In this paper we report for the first time a photonic bandgap fiber-based spectrometer that 
consists of a photonic cystal fiber bimdle and a black and white CCD camera. 

2. Characteristics of the subcomponents: Bragg fibers, fiber bundle, and a CCD camera 

2.1 Photonic bandgap Bragg fibers 



Photonic bandgap fiber is a key element of our spectroscopic system. The photonic 
bandgap Bragg fibers used in our research are home-made and have been reported previously 
in [10-12]. An individual Bragg fiber features a large 300um-700um diameter core made of a 
PMMA plastic. The core region is surrounded with a periodic multilayer reflector featuring 
-100 submicron-thick layers of low and high refractive index PMMA/Polystyrene plastics 
with refractive indices 1.49/1.59. Such a multilayer (Bragg reflector) is responsible for the 
appearance of the spectrally narrow transmission band (reflector bandgap) within which the 
light is strongly confined inside the fiber core. For the wavelengths outside of the reflector 
bandgap the light penetrates deeply into the multilayer region exhibiting strong propagation 
loss due to scattering on the imperfections inside the multilayer region. A typical fiber loss 
within the reflector bandgap region is -lOdB/m and is mostly determined by the bulk 
absorption loss of a low-purity PMMA plastic. Outside of a bandgap region scattering loss 
dominates resulting in >60dB/m propagation loss. In our experiments we used 30cm-long 
fibers so that the loss of guided light was below 3dB, while the loss of non-guided hght 
was >20dB. Numerical aperture of all the Bragg fibers was in the range of 0.17-0.22. 

To construct a fiber bundle spectrometer we chose 100 Bragg fibers with complimentary 
and partially overlapping bandgaps as shown in Fig. 1(a). All the fibers in a bundle were 
drawn from the same preform with the only difference between them being the final diameter. 
The smaller diameter fibers feature bangaps shifted towards the blue part of a spectrum. In 
Fig. 1(b) we also show distribution of the fiber bandgap width as a function of the fiber 
bandgap central wavelength. We note that the fiber bandgap size increases with its center 
Irequency. Thus, in the blue part of a spectrum the fiber bandgap width is ~60nm, while 
increasing to ~180nm in the red part of a spectrum. 
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Figure 1 . (a) Transmission spectra of 5 typical Bragg fibers used in the fiber bundle, (b) Distribution of the fiber 
bandgap width as a flinction of the fiber bandgap center wavelength for all the Bragg fibers in a bundle. In the inset: 
photo of a Bragg fiber crossection showing a solid core surrounded by a periodic multilayer reflector. When 
launching white light into the Bragg fiber the non-guided colors are strongly iiradiated with only a single color 

reaching the fiber end. 



2.2 Photonic bandgap fiber bundle 

The fiber bimdle used in our experiments featured a 5.6mm inner diameter plastic tube 
that hosted 100 Bragg fibers of 30 cm of length. Fibers at the input end of a bimdle were 
rigidly attached with epoxy to each other and to the confining tube, and the whole assembly 
was then polished using optical films of various granularity. On the other end of a bundle 100 
fibers were inserted into a custom-made block featuring holes of diameter 0.8mm placed in a 
periodic IQxlO square array, see Fig. 2. The block, the tube, and all the fibers were bonded 
together with epoxy, and the output end was then polished. 

Principle of operation of a photonic bandgap fiber bundle spectrometer can be clearly 
imderstood from Fig. 1 and Fig. 2. In principle, if all the fiber bandgaps were to be spectrally 
narrow and non-overlapping, then at the end of a fiber bundle the relative intensities of hght 
coming out of the individual fibers would be unambiguously related to the corresponding 
spectral components of an incoming light. In practice, individual fiber bandgaps are always 
overlapping, and in our particular implementation the bandgaps are quite broad. Therefore, to 
reconstruct the intensity of an incoming hght from the intensities of light coming out of the 
individual fibers of a bundle we have to use a certain deconvolution algorithm. Note that 
spectrometer resolution is not directly limited by the width of fiber bandgaps. As we will 
demonstrate later in the paper, even with all the individual bandgaps as wide as 60nm one 
could, in principle, reconstruct the peaks as narrow as 5-1 Onm. However, to achieve such a 
fine resolution, management of the experimental becomes critical. 




Figure 2. Fiber bundle spectrometer. Top part: schematic of the spectrometer. Light from the illuminant is launched 
into the fiber bundle. Lower part: when the broadband light is launched into the fiber bundle, the output is a mosaic 
of colors selected by the individual Bragg fibers. The image is then recorded by the black and white CCD array. 



2.3 Sensitivity and linear response of a CCD sensor 

The output of the fiber bundle is a 10x10 matrix of coloured fibers which is then projected 
with a beam compressor onto a black and white CCD array (see Fig. 2). The sensitivity of a 
CCD sensor would directly affect the measuring range, the signal-to-noise ratio, as well as the 
algorithm that we use to reconstruct the spectrum of an illuminant. The CCD sensor in our 
experiment is an Opteon Depict 1, E and S Series BIA 652x494 Black & White CCD. The 
normalized sensitivity for different wavelengths is displayed in Fig. 3(a) from which it is 
clear that the sensor covers well all the visible and a part of a near-IR spectral range. 
Sensitivity curve, however, is a highly variable function of wavelength, which should be 
taken into account in the re -construction algorithm. Fig. 3(a) was obtained by launching the 
light form a supercontinuum (SC) source into a monochromator whose output would produce 
a 2nm-wide (FWHM) peak at a desired center wavelength. Light from such a tuneable source 
would then be launched into a commercial multimode fiber of ~lmm diameter. The output 
power from such a fiber would then be measured by both the CCD array and a calibrated 
power meter; the sensitivity curve is then achieved by dividing one by the other 



As we will see later in the presentation of our paper, one of the key requirements to a 
CCD array is linearity of its monochromatic response as a function of the intensity of incident 
light. The near-linear response of a CCD sensor was confirmed at several wavelengths of 
interest; a typical CCD response (at ?L=560nm) is presented in Fig. 3(b). To obtain Fig. 3(b), 
two polarisers were placed between a CCD array and the output of a commercial fiber in a 
setup described in the previous paragraph. By varying the angle between the two polarizes we 
also varied the intensity of light coming onto a CCD array. A reference measurement is 
performed with a calibrated power meter for the same angles between polarisers; the results 
are compared to establish a near-linear sensor response. 
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Figure 3. (a). Normalized spectral response of a CCD array, (b) A typical monochromatic near-linear response 
of a CCD array (?i=560 nm). (c) An image of fiber bundle taken by a CCD array. 

3. Calibration of tlie fiber bundle spectrometer and spectrum reconstruction algorithm 

In its operation mode the fiber bundle spectrometer is illuminated with a test light, which 
is then spectrally and spatially decomposed by the fiber bundle see (see Fig. 4). The test 
spectra used in our experiments were either broadband (halogen lamp source), or narrowband 
(tuneable monochromator source). Image of the fiber bundle output end was then projected 
onto a Black and White CCD array using a 2:1 beam compressor. Note that a CCD array does 
not have to be a 2D camera, instead one can use a much economical linear array with the 
number of detectors matching the number of Bragg fibers in a bundle. Output of each fiber 
carries information about the intensity of a certain fraction of an illuminant spectrum as 
filtering by the individual Bragg fibers. It is reasonable to assume then that if the Bragg fibers 
of a bundle feature complimentary bandgaps that together cover the spectral range of interest, 
then the spectrum of an illuminant could be reconstructed. 



Fig.4. Setup for the spectrometer calibration measurement. Indicated from left to right are the tuneable narrowband 
monochromator-based source, a difluser to create a uniform illuminance of the fiber bundle input end, photonic 
bandgap fiber bundle, a beam compressor made of the two biconvex lens, and a CCD sensor. 

3. 1 Transmission matrix method 

To interpret B&W CCD image, we use a so-called luminance adaptation model of the 

fiber spectrometer. Particularly, for a fixed exposure time of a CCD array, total intensity 

registered by a CCD sensor from the output of the n* fiber can be presented as: 

c„ = \';'^ mAFsi)s{X)osm^ (1), 

where /(/I) is the illuminant spectral flux density at the fiber-bundle input end, ^„ is 
the effective area of the n"' fiber, F^i^X) is the transmission function of the n* fiber, S(^X) 

is the spectral sensitivity of a CCD array, 0„ (/I) is a fiber-position dependent transmission 

function of various optics (diffiiser, beam compressor, and a CCD objective). Measuring such 
transmission functions individually is a daunting task. Instead, we use a calibration procedure 
that measures a compounded transmission matrix of our spectrometer defined as 

7^, (/I) = A^F^ (/I) . The model equation (1) could then be now rewritten as: 



A.- +A/i 



C„ = r^IiX)T„iX)dA= X TS^) \ mdA = J^I.Tl (2). 



!=0 



Experimentally, an integral in equation (2) is rather a sum over the intensity contributions 
from small spectral "bins" of size A/l . Equation (2) presented in the mafrix fr)rm becomes: 
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where vector (C^^^j^^j represents the intensities of light coming out of the individual 
fibers as measured by the CCD sensor, is a spectrometer transmission mafrix, and 

^^^Nxi ^ discretized spectrum of the illuminant. From (3) it follows that a discretized 

spectrum of the illuminant can be reconsfructed from the corresponding CCD image (C 
vector) by inverting the transmission matrix of a specfrometer. 

3.2. Calibration measurement, building a transmission matrix 

To construct transmission matrix experimentally we note that if the illuminant spectrum is 
monochromatic ( I^ = 0, V/ in (3)), then the measured C'^ vector is proportional to the 

column of the fransmission mafrix: 

CHl:100) = r(l:100,iJ-/. (4). 

To construct the transmission matrix experimentally we use a tuneable 
monochromator-based narrow-band (2nm FWHM) source to generated "monochromatic" 
specfra (see Fig. 4). Particularly, we vary the source center wavelength in 2nm increments, 
thus effectively subdividing the 400nm-840nm spectral interval under consideration into 
N = 22l equivalent 2nm-wide bins. For every new position of the source center 
wavelength we acquire a C vector using a CCD array and consider it as the next column of 
the specfrometer transmission matrix. Finally, to finish calibration we measure the 

wavelength -dependent intensities of a tuneable source by placing a calibrated power 

meter directly at the output of a monochromator By dividing every C vector by the 

corresponding / value the transmission matrix is constructed. 

In our experiments an S&Y halogen lamp source was used with a Newport Oriel l/8m 
monochromator to build a narrow-band tuneable source of 2.0nm FWHM (see Fig. 4). The 
light beam from a source was then directed onto a diffuser placed right before the fiber 
bundle to guarantee a uniform illumination of its input end. At the output end of a fiber 



bundle a 2:1 beam compressor was used to image the fiber bundle output facet onto an 8-bit 
B&W CCD sensor array. Individual images were then interpreted using a Matlab image 
processing suite to construct the corresponding C vectors. Wavelength dependent intensity 
of a tuneable source was measured using a calibrated Newport 841-PE power meter. 

3.3 Spectral reconstruction algorithm 

As it was noted earlier, a discretized spectrum of the illuminant can be reconstructed from 
the corresponding CCD image ( C vector) by inverting the transmission matrix of a 
spectrometer in (3). However, an immediate problem that one encounters when trying to 
invert the transmission matrix is that the matrix is non-square. Even if the number of spectral 
bins is chosen to match the number of fibers in the bundle, thus resulting in a square 
transmission matrix, one finds that such a matrix is ill-conditioned. One, therefore, has to 
resort to an approximate inverse of a transmission matrix. To find a pseudo-inverse of a 
transmission matrix in equation (3) we employ a singular value decomposition (SVD) 
algorithm. Particularly, form the linear algebra we know that any (lOOxN) matrix (suppose 
that N>100) can be presented in the form: 

(t) =(u) (s) (v^) 

V /lOOxA- V /lOOxlOO V /IOOxjV V /A/xA' 

U^U = 1 ; V^V = \ ; S = diag(a„a^,a„.. .a,^ ) (4), 

o-j >o-2 >o-3 >...>o-ioo >0 

where matrixes U and V are unitary, and matrix 5* is a diagonal matrix of the real 
positive singular values. For the ill-conditioned matrices, most of the singular values are 

small and can be taken as zero. By limiting the number of non-zero singular values to A'^^ , 

equation (4) can be rewritten as: 

("^)lOOxiV ~ (^)lOOxA'^ i^^N^xN^ xAf 



cr,,cr2,...,o-^ 



(5), 



and a corresponding pseudo-inverse of the transmission matrix is then found as: 

xlOO 



(-^)AfxlOO (^)iVxiV, i^)N„xN„ )w xlOO 



5. Experimental results 

To test our fiber bundle spectrometer we consider resolving a set of the 25nm-wide peaks 
centered at 450nm, 500nm, 550nm, 600nm, and 700nm. Such peaks were created using the 
same monochromator-based tuneable source as used for the spectrometer calibration, 
however adjusted to have a 25nm bandwidth of the outgoing light. In Fig. 5 we demonstrate 
in black dashed curves the test spectra of light beams at the input of a fiber bundle (as 
resolved by another Oriel monocromator), while in thick red lines we show the corresponding 



reconstructed spectra using the transmission matrix inversion algorithm. In the figures we 
also indicate the optimal number of non-zero singular values used in the matrix inversion 
procedure (6). The choice of the number of non-zero singular values used in the spectrum 
reconstruction algorithm effects strongly the quality of a reconstructed spectrum. Note, in 
particular, that although the spectral intensity function should be strictly non-negative, the 
reconstructed spectral intensity can take negative values. Therefore, spectrum reconstruction 
error can be defined as the ratio of the most negative value of the reconstructed spectral 
intensity to its most positive value: 

Error{N^) = -xmn{I~) I (7), 



and the optimal A^^ to be used in (6) is the one that minimises the error (7). 
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Figure 5. Spectra reconstruction using photonic bandgap fiber bundle-based spectrometer In black dashed curves are 
the test spectra of the incoming light as created by the tuneable monochromator-based source, hi red solid curves are 

the corresponding spectra reconstructed by the fiber bundle spectrometer. No indicates the optimal number of 
singular values used in the transmission matrix inversion algorithm. Gray areas indicate error levels. 

From Fig. 5 we conclude that the test peak position (the value of the peak center 
wavelength) can be reconstructed within several percent of its true value in the whole spectral 
rage covered by the spectrometer. Moreover, we note that although the widths of the 
individual Bragg fiber bandgaps are quite large (>60nm) the fiber bundle-based spectroscopic 
system can resolve peaks of much smaller spectral width. Indeed, in Fig. 5, the widths of the 
reconstructed spectra are in the 30-50nm range with the exception of a peak centered at 
600nm. In the next section we present theoretical estimates for the spectral resolution limit of 



our fiber bundle spectrometer. We mention in passing that we have also conducted 
measurements using 5nm, lOnm, 15nm, and 20nm-wide peaks. In all the experiments we saw 
that the center wavelength of even the narrowest peak can be reconstructed with high 
precision, however the reconstructed peak width always stayed larger than 30nm. 



6. Limits of spectral resolution for the fiber bundle-based spectrometer 

To understand the resolution limit of our spectrometer we first study the effect of the 

choice of the number of non-zero singular values A^^ used in the inversion algorithm (6). At 

first we assume that there is no noise in the system. As a test spectrum /(/I) we consider a 
20nm-wide peak centered at the wavelength of 520nm (dashed black curve in Fig 6(a)). 
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Figure 6. Properties of the reconstructed spectra as a fiinction of the number of singular values used in the inversion 
algorithm. No noise is present in the system, (a) Dependence of the spectral shape of a reconstructed peak on A'^ . (b) 



Width of a reconstructed peak as a fUnction of A'^ . (c) Reconstruction error as a function of N . 

We then multiply the discritized test spectrum by the transmission matrix of our 
spectrometer (measured experimentally), and then reconstruct the test spectrum by using a 

transmission matrix pseudo-inverse with A'^^ singular values: 

(^)lOOxl ~ ("^)lOOxA' (^)iVxl 



Thus reconstructed spectrum is then compared to the original test spectrum. In Fig. 6(a) 
we present several reconstructed spectra for the different number of singular values used in 



the inversion algorithm. If all the 100 singular values are used (red thin curve in Fig. 6(a)), 
then the 20nm-wide peak is very well reconstructed featuring less than 5% intensity 
difference with the reconstructed spectrum. The error manifests itself in the form of ripples at 

the peak tails. When smaller number of singular values are used (for example A^^. =10, 

blue solid cure in Fig. 6(a)), center wavelength of the reconstructed peak still coincides very 
well with that of a test peak, however, reconstructed peak width is larger than that of a narrow 
test peak. This result is easy to understand by remembering that expansion basis set used in 
the reconstruction algorithm is formed by the relatively broad (>60nm see Fig. 1) 
transmission spectra of the Bragg fibers, therefore one needs a large number of such basis 
functions (singular values) to reconstruct a narrow spectral feature. In Fig. 6(b) the FWHM of 
the reconstructed peak is shown as a function of the number of singular values used in the 
matrix inversion. Note that more than 40 singular values have to be used to achieve the width 
of a reconstructed peak to be close to 20nm width of a test peak. When only 10 singular 
values are used, the width of a reconstructed peak increases to 40nm. Additionally, the peak 
reconstruction error as defined by (7) increases to 20% when a small number of singular 

values A^^ < 20 is used. From this we conclude that with no noise present in the system, 

the peak reconstruction error and the quality of a reconstructed spectrum improves when a 
larger number of singular values is used. 

We now study the effect of experimental noise on the quality of a reconstructed spectrum. 
There are several sources of noise in our spectrometer. One is the discrete intensity resolution 
of a CCD amy; with the 8-bit encoding the fundamental noise is Moreover, because of 
the relatively large size of a fiber bundle (~6mm) it is difficult to ensure the same 
illumination conditions of the bundle input facet for all the sources used in the experiments. 

1.4 

I ''^ 

S 1 

=! 
O 
U 

'^0,8 
So,6 

Q. 

£ 0.4 
o 
o 

t 0.2 






(III) fl 




6-0.0S,N„=15 

Test specdTim 









6=0.1 








f 


K . 


aio.o5 

.„ .■' 8=0.005 


(ivf (iH)^ 














'J 


SK>;Nfl=100 

Test spectrum 







600 700 
X(niTi} 



40 60 



600 700 



Figure 7. Effect of noise on tlie quality of reconstruction. Examples of the reconstructed spectra for several particular 
realisations of noise with amplitudes: (i) ^ = , (ii) S = 0.005 , (iii) S = 0.05 , (iv) S = 0.1 . 

To model experimental noise numerically, we first multiply the test spectrum by the 
transmission matrix of our spectrometer, then add a uniformly distributed random noise of the 



relative intensity S to the C vector (B&W image), and finally reconstruct the test 
spectrum by using a transmission matrix pseudo-inverse with A^^ singular values: 



100x1 V ^lOOxAT V JNxI 

(^""Li=(CLxr(l + ^-('7)iooJ ; rj- random ^[-0.5,0.5] (9). 

(Treconstr.\ Tl^^ fiS*^"^ (U^\ f ^ noise \ 

With these definitions, we first revisit the case when all the 100 singular values are used 
for spectrum reconstruction. As noted above, in the absence of noise ( = in Fig. 7(i)) the 
20nm peak can be well reconstructed with only ~5% error. However, addition of even a small 
amount of noise {d = 0.005 in Fig. 7(ii)) results in a highly noisy reconstructed image with 
errors as large as 50%. Note that by using a smaller number of singular values 

A/^(= [20,40] (see the error plot in the middle of Fig. 7) the error of the peak 

reconstruction can be greatly reduced down to 10%, while still allowing a fair estimate of the 
peak width ~30nm. When the noise level is increased further, to avoid large reconstruction 

errors one has to use a relatively small number of singular values < 20 (see Figs. 

7(iii,iv)). Notably, reconstruction error saturates at 20% even for large noise levels, however, 
the reconstructed peak width stays always larger than 30nm due to small number of singular 
values used. 

Finally, Fig. 8 presents statistical averages of various parameters (thick sohd lines), as 
well as their statistical deviation irom the average (thin dashed curves) as a function of the 
noise level. To construct these curves, for every value of the noise amplitude S we, first, 
generate 200 realizations of noise. Then, for every realisation of noise we plot the error of 
peak reconstruction as a function of the number of singular values (similar to Fig. 7). Lastly, 

the smallest error, the corresponding and the width of a reconstructed peak are recorded 

and statistical averages are performed. In Fig. 8(a) we present the average peak reconstruction 
error as a function of the noise amplitude. Note that although the error grows with the noise 
amplitude, it, nevertheless, saturates at -20%) even for large noise levels. In the inset of Fig. 
8(a) we see that at very small noise levels {5 < 0.02 ) to obtain the smallest reconstruction 
error one can use almost all 100 singuleir values, while at higher noise levels less than 20 
singular values should be used. As a consequence, for larger noise amplitudes the width of a 
reconstructed peak could become considerably larger than that of a test peak (20nm). 
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Figure 8. Effect of noise on tlie reconstruction algorithm, (a) Average reconstruction error and its statistical deviation 
as a function of the noise amplitude. Inset: optimal number of singular values needed to minimise the reconstruction 
eiTor (b) Average width of a reconstructed peak and its statistical deviation as a function of the noise level. 

7. Conclusion 

We have demonstrated experimentally a novel all-fiber spectrometer comprising a 
photonic bandgap fiber bundle and a black and white CCD camera. 100 Bragg fibers with 
complimentary bandgaps covering a 400nm-840nm spectral range were placed together to 
make a ~6mm in diameter, 30cm-long fiber bundle. The spectrometer operates by launching 
the light into a fiber bundle and then recording the color-separated image at the fiber bundle 
output facet using a CCD camera. Inversion algorithm was developed to reconstruct the test 
spectrum using black and white intensity images. Among the clear advantages of our 
spectrometer is the lack of moving parts, a near instantaneous and parallel acquisition of all 
the spectral components, compactness, high degree of integration and simplicity of operation. 

Photonic bandgap fiber bundle spectrometer was demonstrated experimentally to resolve 
consistently and with high precision the center wavelengths (positions) of the relatively 
narrow (5nm-25nm) spectral peaks. If the peak width is much narrower than the bandwidth of 
Bragg fibers (>60nm) making the fiber bundle, the peak width can still be resolved correctly, 
however, the reconstruction algorithm is sensitive to experimental noise. Experimental 
resolution limit of the fiber bundle spectrometer presented in this paper is found to be ~30nm. 
Theoretical analysis predicts that resolution limit of the existing setup could be greatly 
improved via careful management of the experimental errors. 
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